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ABSTRACT 

We present a detailed examination of the brown dwarf multiples 2MASS J08503593+1057156 and 
2MASS J17281150+3948593, both suspected of harboring components that straddle the L dwarf/T 
dwarf transition. Resolved photometry from Hubble Space TeZescope/NICMOS show opposite trends 
in the relative colors of the components, with the secondary of 2MASS J0850+1057 being redder than 
its primary, while that of 2MASS J1728+3948 is bluer. We determine near-infrared component types 
by matching combined-light, near-infrared spectral data to binary templates, with component spectra 
scaled to resolved NICMOS and K p photometry. Combinations of L7 + L6 for 2MASS J0850+1057 and 
L5 + L6.5 for 2MASS J1728+3948 are inferred. Remarkably, the primary of 2MASS J0850+1057 ap- 
pears to have a later-type classification compared to its secondary, despite being 0.8-1.2 mag brighter 
in the near-infrared, while the primary of 2MASS J1728+3948 is unusually early for its combined-light 
optical classification. Comparison to absolute magnitude/spectral type trends also distinguishes these 
components, with 2MASS J0850+1057A being wl mag brighter and 2MASS J1728+3948A «0.5 mag 
fainter than equivalently-classified field counterparts. We deduce that thick condensate clouds are 
likely responsible for the unusual properties of 2MASS J1728+3948A, while 2MASS J0850+1057A is 
either an inflated young brown dwarf or a tight unresolved binary, making it potentially part of a 
wide, low-mass, hierarchical quintuple system. 

Subject headings: binaries: visual stars: individual (2MASS J08503593+1057156, 

2MASS J17281150+3948593) — stars: low mass, brown dwarfs 



1. INTRODUCTION 

The L dwarfs and the T dwarfs are the two low- 
est luminosity spectral classes of very low-mass stars 
and brown dwarfs known today, encompassing effec- 
tive temp eratures (T e ff) f r om ~2200 K down to 
~600 K (IGolimowski et all 120041 : iBurningham et al. 
120081: iStephens et al.l 120091: see also iKirkpatrickl I200E 
and references therein). Spectroscopic studies of these 
sources reveal atmospheres that are remarkably di- 
verse. L dwarfs having abundant molecular gas species 
and clouds of condensates in their photospheres, while 
T dwarfs have relatively cloud-free photospheres and 
more complex molecular gas species including CH4, 
NH3 and CO9 (lOppenheimer et ai l 119951: iRo ellig et al.l 
120041 : lYamamura et al.l 120101 ). Condensate cloud prop- 
erties are believed responsible for near-infrared spec- 
tral and c olor variation s among equivalently-c l assified 
L dwarfs (|Knapp et al.l I2004t IBurgasser etal] I2008bt 
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ILooper et al.l I2008bf) . temporal variability in late-type 
dwarf s (IBailer- Jones fc Mundtl [200lt IGelino et all 12001 
iKoenl [20031 : lArtigau et al.l 12009ft and dramatic changes 
in the spectral energy distributions between the L dwarf 
and T dwarf classes (iLeggett et all 120001: lAllard et al.l 
[20011 lAckerman fc Marlevi 120011: IBurrows et al.l 12006). 
Non-eq uilibrium chemistry and atmospheric circu- 
lation (IBurrows et al.l 120001: ILodders fc Feglevl 120021: 
iSaumon et al.ll2006|) coupled with the complex processes 
of condensat e grain formation, growth, circulatio n and 
evaporation ([Helling fc Woitkdl2006l: [Helling et al.ll2008T) 
makes the formation, evolution and influence of conden- 
sate clouds in low-temperature atmospheres among the 
outstanding problems in brown dwarf astrophysics today. 

The disappearance of photospheric condensates at the 
transition between the L dwarf and T dwarf classes is 
one particularly interesting aspect of brown dwarf clouds. 
This transition occurs over a narrow range of tempera- 
tures (A T e ff ss 200 K) and luminosities (AlogLf, ;/L(^ 
0.2 dex; [K irkpatri ck et al.ll2000t IGolimowski et al.H2004l 
lBurgasserll2007afl . and is accompanied by a temporary 
brightening at 1 fj,m (the "J-band bump " ; iDahn et all 
[200alTinnev et al.| [2003l: iVrba et al.ll2004ri an d enhanced 
rates of multiplicity (jBurgasser et al.l 120063 ). These ef- 
fects appear to arise from the rapid depletion of conden- 
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sate clouds at the L dw arf/ T dwarf transition ([Liu et al.l 
120061 : [Burgasser 2007a), although the driving mechanism 
for that deple tion is inadequately explained by current 
cloud models (iBurrows et all 120061 iCushing et al.1 120081 : 
iSaumon fe Marlevl 120081 1 Whether arising from the 
fragmentation of cloud stru ctures (jAckerman fe Marlevl 
120011 : [Burgass er et al.l [2002T) . a sudden increase in sedi- 
ment ation efficiency (jKnapp et al.l [2004: Ste phens et al.l 
2009) or some other process remains unclear, but cloud 
depletion at the L dwarf/T dwarf transition has im- 
portant implications on cloud evolution in other low- 
tempera ture atmosph eres, such as those of extrasolar 
planets ([Fortnevi 120051 ) . 

Coeval and cospatial multiples are important labo- 
ratories for studying this transition, eliminating de- 
pendencies on age, composition and distance. So- 
called "flux-reversal" binaries, whose components strad- 
dle the L dwarf/T dwarf transition, have verified the 
1 fim brightening as an intrinsic aspect of brown 
dwarf evolution (IGizis et al.l l200l iCruz et al.l 12004 
Burgasser et al.l l2006ct iLiu et al.l 120061 : lLooper et al.l 
2008aj) . Yet resolved spectroscopy of the compo- 
nents of these systems has been difficult to obtain, 
due largely to the close separa tions typica l of brown 
dwarf multiples (a < 20 AU; lAllenl (20071) . In this 
article, we present a comprehensive analysis of two 
binaries whose components are suspected to straddle 
this transition: 2MASS J08503593+1057156 (hereafter 
2MASS J0850+1057) and 2MASS J17281150+3948593 
(hereafter 2MASS J1728+3948). Both were resolved 
with the Hubble Space Tele scope (HST) Wide Field Plan- 
etary Camera 2 (WFPC2; IReid et al.ll200lB: IGizis et al.l 
|2003|) and noted as potential L dwarf/T dwarf tran- 
sition binaries based on their photometric properties 
and late (combined-light) systemic spectral types of L6 
and L7, respectively. Here we combine resolved near- 
infrared photometry obtained with the HST/Near In- 
frared Camera and Multi-Object Spectrometer (NIC- 
MOS) and ground-based laser guide star adaptive optics 
(LGSAO) imaging with combined-light, near-infrared 
spectroscopy to infer component spectral types, colors 
and absolute magnitudes. Observations are described 
in Section 2. In Section 3 we describe point spread 
function (PSF) fitting of the HST data that yield rela- 
tive fluxes for each binary in five spectral bands span- 
ning 1.0-1.8 /im. In Section 4 we describe our spec- 
tral fitting procedure and determine component classi- 
fications, colors and absolute magnitudes. In Section 5 
we compare these measures to current absolute magni- 
tude/spectral type and absolute magnitude/color trends. 
In Section 6 we discuss our results in the context of cloud 
evolution at the L dwarf/T dwarf transition, and explore 
the possibility of higher-order multiplicity in the case of 
2MASS J0850+1057. Results are summarized in Sec- 
tion 7. 

2. OBSERVATIONS 
2.1. Targets 

The empirical properties of the two binaries examined 
here are summarized in Table ffl 2MASS J0850+1057 
was initially identified in the Two Mi cron All Sky Sur- 
yey_(2MASS: lSkrutskie et al.ll2006l) bv lKirkpatrick et al.l 
(1999), and selected to be the prototype for the L6 



spectral subclass. Its optical spectrum exhibits 6708 A 
Li I absorp tion, indicating component masses below 
~0.06 Mq (|Rebolo et al.l 119921: iMagazzu et all [19931) . 
This sourc e was r esolved as a 0'.'16 binary system by 
IReid et al.l (|2001bD with HST/WFPC2, and has been 
subsequently confirmed as a common proper motion 
pair from multi-epoch HST and LGSA O observations 
([Bouv et al.ll2008HKonopackv et al.ll20Tol) . The large dif- 
ference in component brightne sses in the WFPC 2 F814W 
band (AF814W = 1.47±0.09: lBouv et al.ll2003T) has sug- 
gested a late-type L or early-type T dwarf secondary, 
although this component has been a persistent outlier 
in color magn itude diagrams. As trom etric parallax mea - 
surements by IDahn et all (|2002j) and IVrba et al.l (|2004l ) 
find discrepant distan ces of 25.6±2.3 pc and 38±6 pc, dif- 
fering by nearly 2a. iFahertv et al.l (|2010f ) identified an 
unrelated b a ckgro und source that may have skewed the 
IDahn et al.l (|2002f ) astrometric measurements, and re- 
port a preliminary parallax that is intermedia t e betw een 
these two values. We adopt t he IVrba et al.l (120041) fo r 
this study, as was al so used bv lKonopackv et al.l (|2010l ). 
IFahertv et al.l (|2010T ) also identified a widely-separated, 
common proper motion and common distance compan- 
ion to 2MASS J0850+1057, the M5 + M6 binary NLTT 
20346AB at a projected separation of 4.'1 (7700 AU). Ha 
and X-ray emission in these M dwarfs, coupled with Li I 
absorption in 2MASS J0850+1057, indicate a relatively 
young age of 0.25-1.5 Gyr for the combined system, and 
relatively low masses for the 2MASS J0850+1057 com- 
ponents: 0.04±0.02 M w and 0.03±0.01 M^ based on 
evolutionary models (|Fahertv et al.l [2010). These mass 
est imates are consist e nt wit h weak empirical constraints 
bv iKonopackv et all (|2010l ). M tota; = 0.2±0.2 M , as 
inferred from incomplete astrometric monitoring of its 
orbit. 

2MASS J1728+3948 w as also identified in 2MASS by 
iKirkpatrick et ail ([2000D and classified L7 based on its 
optical spectrum. This source shows no indication of 
Li I absorption and hence its pri mary is infe r red to have 
a mass gre a ter th an 0.06 M Q . IGizis et"aT1 ([20031 ) and 
IBouv et al.l ([20031 ) identified this system as a 0'.'13 bi- 
nary based on HST/WFPC2 observations. The former 
study found the "secondary" to be fainter in the F814W 
band and brighter in the F1042M band, the first reported 
example of a L/T flux reversal binary. Multi-epoch HST 
and LGSAO observations rep orted in lBouv et al.l (2008) 
and IKonopackv et all (|2010l ) have verified the common 
proper motion of the components, with the latter study 
finding an orbital period of 31±12 yr, semimajor axis of 
5.3±0.8 AU (ba sed on the paralla x distance measurement 
of 24.1±1.9 pc; IVrba et al.l [2004) and total system mass 
of 0.15t£o4 M . However, like 2MASS J0850+1057, the 
orbit of this system has not been fully sampled. As- 
suming a primary temperature of «1450 K based on 
the c o mbined light optical spe ctral type ([Stephens et al.l 
120091: IKonopackv et all I2010D . the absence of Li I in- 
dicates a system age of at l east 1.5 Gyr, accordi ng to 
the evolutionary models of iBurrows et al.l ([1997D and 
ISaumon fc Marlevl ((20081) . 

2.2. IRTF/SpeX Observations 

Low resolution near-infrared spectral data for 
2MASS J0850+1057 and 2MASS J1728+3948 were 
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TABLE 1 

Properties of 2 MASS J08503593+1057156 and 2MASS J17281150+3948593 



Parameter 2MASS J0850+1057 2MASS J1728+3948 Reference 



Optical Spectral Type a 


L6 


L7 


1,2 


NIR Spectral Type a > b 


L7±2 


L6±l 


3 


MKO J 


16.20±0.03 


15.90±0.08 c 


3,4 


MKO J — K 


1.85±0.04 


2.01±0.09 c 


3,1 


Distance (pc) 


38±6 d 


24.1±1.9 


5 


Vtan (km s _1 ) 


26.6±4.5 


5.1±0.9 


5 


Separation (mas) e 


157±3 


131±3 


6 


(AU) 


6.0±0.9 


3.2±0.3 


5.6 


Position Angle (°) e 


114.7±0.3 


27.6±1.2 


6 


AF814VV 


1.47±0.09 


0.37±0.04 f 


6.7 


AF1042A/ 




-0.25±0.14 


7 


AK p s 


0.78±0.07 


0.63±0.03 


8 


Li I? 


Yes 


No 


1.2 


Age (Gyr) h 


0.25-1.5 


>l-5 


3,9 


Combined Mass (M )' 


0.2±0.2 


15+ 025 
U - 1O -0.04 


8 



Referen ces. — (1) |Kir kpatrick ct al. (1999): (2) Kirkpatrick ot al. (20 0(3); (3) This 
paper; (4) ISkrutskie et al.l (1200611 : (51 IVrba et all 1200411 1; f 61 IBouv et al.l 0200311 ; (7) 
IGizis et all 120031 : (81 IKonopackv et alj (1201011 ; (9) IFahertv et al.l (1201011 . 
a Classification of combined light spectra. 

Based on spectral indices measured from SpeX spectroscopy; see Section 2.2 and Ta- 
bic [2] 

c Synthesized from 2MASS photometry and spectrophotometry filter corrections com- 
puted from the spectrum shown in Figure [T1 

Dahn ct al. (2002) and Fahcrtv ct al. (2010J) report parallactic distances of 25.6±2.3 pc 
and 29±7 pc for 2MASS J0850+1057, both closer than the IVrba et al.l ( 1200411 measure- 
ment. The Fahcrty et al. study propose that the differences may arises from a contaminant 
background source skewing ccntcr-of-light measurements. We adopt the Vrba et al. dis- 
tance for consistency, and note that an uncertainty-weighted mean of this measurement 
and that of Fahcrty et al., 36±5 pc, is fully consistent with the adopted value. 
° At epoch 2000 February 1 (UT) for 2MASS J0850+1057 and 2000 August 12 (UT) from 
2 MASS J1728+3948 . 

f IBouv et"aT1 J2003I1 report AF814W = 0.66±0.11 for this system based on the same 
dataset, a 2.5a" difference. 

s Uncertainty-weighted averages of multi-epoch flux ratio measurements from 
Konopacky ct al. (2010), using both statistical and PSF-matching systematic uncertain- 
ties. 

Based on presence/absence of Li I, and for 2MASS J0850+1057 activity diagnostics in 
its co-moving companion NLTT 20346AB dFahertv et al.H2010H . 
1 Based on partial coverage of astromctric orbits. 



obtained with the 3m NAS A Infrared Telescop e Facility 
(IRTF) SpeX spectrograph (jRavner et al.|[2003h on 2008 
January 8 (UT) and 2006 August 21 (UT), respec- 
tively. Conditions on both nights were mostly clear, 
with light cirrus during the 2006 August observations; 
seeing was ^0'.'6 at J-band in both observations. We 
employed the prism-dispersed mode of SpeX which 
provides 0.75-2.5 /mi continuous spectroscopy at a 
resolution A/AA « 120 for the 0V5 slit, and dispersion of 
20-30 A pixeU 1 . The slit was aligned to the parallactic 
angle in all observations. 2MASS J0850+1057 was 
observed at an airmass of 1.13, with four exposures of 
150 s each obtained in an ABBA dither pattern along 
the slit. 2MASS J 1728+3948 was observed in a similar 
manner, at an airmass of 1.06 and with four exposures 
of 180 s each. After each target observation, AO V 
stars HD 74721 (V = 8.71) and HD 164899 (V = 7.91) 
were observed with identical instrument settings and 
at a similar airmass for telluric absorption and flux 
calibration. These were followed by internal flat field 
and Ar arc lamps for pixel response and wavelength cali- 
bration. Da ta were reduced using the SpeXtool pa ckage, 
version 3.3 (TVacca et al.|[2u0llCushing et alJl200l using 
standard settings; see iBurgasser fc McElwainl (|2006f ) for 



details. 

The reduced spectra for both sources are shown in Fig- 
ure [U The data exhibit classic signatures of L dwarf 
near-infrared spectra, including a steep 0.8-1.2 /xm spec- 
tral slope; FeH absorption at 1.0 /mi; unresolved K I 
and Na I doublets at 1.13, 1.17 and 1.22 /jm; deep H 2 
absorption bands cent ered at 1.4 and 1.9 ^m; and CO 
absor ption at 2.3 urn (iReid et al.1 l2001at iMcLean et al.l 
120031 : iCushing et al.ll2005f ). Their overall near-infrared 
spectral energy distributions are fairly red, consistent 
with their J — K s colors and indicative of cloud opac- 
ity capping the 1.25 /mi and 1.65 /jm spectral peaks 
(|Ackerman & Marlevj[200ll ). There is no obvious indica- 
tion of CH4 absorption at either 1.6 /mi or 2.2 /mi, com- 
monly seen in th e combined light spectra of L dwa rf/T 
dwarf pairs (e.g. JCruz et al.ll20M lBurgasserll2007bl ). 

We derived near-infrared spectral types for these 
sources using th e H2O-J, H2O-H and CH 4-K spectral in- 
dices defined in IBurgasser et all (l2006bD. and th e index- 
spectral type relations defined in IBurgasser! (|2007al) 
which have a typical scatter of ~1 subtype for L0-L8 
dwarfs. Values are reported in Table [2] The inferred 
near-infrared types are L7±2 for 2MASS J0850+1057 
and L6±l for 2MASS J1728+3948. These types are 
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1.5 2.0 

Wavelength (urn) 

Fig. 1.— Near-infrared spectra of 2MASS J0850+1057 (bottom) 
and 2MASS J1728+3948 (top) obtained with IRTF/Spex. Data 
are normalized at the 1.3 (im spectral peaks, with the spectrum of 
2MASS J1728+3948 vertically offset for clarity (dotted line). Ma- 
jor spectral features characteristic of L dwarf spectra are labeled. 
Also indicated are the wavelength ranges spanned by the F110M, 
F110W, F145M, F160W and F170M HST/NICMOS filters. 



TABLE 2 

Near-Infrared Spectral Index Measurements 





2MASS J0850+1057 


2MASS J1728+3948 


Index 


Value 


SpT 


Value 


SpT 


H 2 0-J 
H2O-H 
CH4-K 


0.673±0.014 
0.652±0.010 
0.987±0.009 


L8.0±0.4 
L8.2±0.4 
L4.8±0.3 


0.712±0.009 
0.723±0.005 
0.955±0.004 


L6.9±0.3 
L5.6±0.2 
L5.9±0.2 


Mean SpT 




L7±2 




L6±l 



Note. — Spectral ind i ces and ind ex-spectral type re lations arc dc- 
fincd in Burgasser ot al. (2006b) and Bui ■gasscr (2007a). respectively. 
Mean values and uncertainties in indices and index types were calcu- 
lated from 1000 realizations of each spectrum randomly varied accord- 
ing its noise spectrum. Overall mean types and their uncertainties arc 
the average and scatter of the index types, rounded off to the nearest 
whole subtype. 

later and earlier than the reported optical types for these 
sources, respectively, albeit consistent within the uncer- 
tainties. The larger uncertainty in the near-infrared spec- 
tral type of 2MASS J0850+1057 is driven by a discrepant 
CH4-K index value, which indicates an L5 spectral type 
as compared to L8 from the H2O indices. This difference 
may be related to the multiplicity of this system, or its 
young age and the corresponding low surface gravity of 
its components. 

2.3. HST/NICMOS Observations 

Both sources were observed in single orbits with the 
HST/NICMOS NIC1 camera as part of program GO- 
9843 (PI Gizis). 2MASS J0850+1057 was observed on 
2003 November 9 (UT) and 2MASS J1728+3948 on 2003 
September 7 (UT). NIC1 is the highest-resolution cam- 
era on NICMOS with pixel scale 0'.'043 and field of view 
ll"xll", providing a well-sampled PSF down to the 
diffraction limit at 1 /tin. Both sources were observed 
through the wide-band filters F110W (A c = 1.025 /im, 



AA = 0.6 /im) and F160W (A c = 1.55 /im, AA = 0.4 /im), 
and the medium-band filters F110M (A c = 1.1 /im, AA 
= 0.2 /im), F145M (A c = 1.45 /im, AA = 0.2 /im), and 
F170M (A c = 1.7 /im, AA = 0.2 /im). As shown in Fig- 
ured! the F110W and F160W filters sample the promi- 
nent J and H flux peaks of late-type L and T dwarf spec- 
tra, while the F110M, F145M and F170M filters sample 
the 1.1 /im H2O, 1.4 /im H2O and 1.6 /im CH4 absorp- 
tion bands, respectively. These filter combinations pro- 
vide discriminating colors for late- type L and T dwarfs 
(e.g- lBurgasser et al.l[2006d iReid et al.ll2006l : Figure [3]). 
All data were acquired in MULTIACCUM mode, with 
two dithered exposures (1'.'5 step) obtained in each fil- 
ter. Total integration times for each source were 256 s in 
F110W and F110M (2MASS J1728+3948 was observed 
for 288 s in F110M), 320 s in F160W, 704 s in F170M 
and 832 s in F145M. Integrations were longest with the 
F145M filter because of the deep H2O absorption present 
at these wavelengths. 

Raw images w ere reduced by s tandard pipeline process- 
ing (CALNICA, Bushousc 1997) using calibration images 
(inflight and model reference files, circa 2002) and pho- 
tometric keywords (post-cyrogenic) current as of Febru- 
ary 2008. CALNICA reduction includes analog-to-digital 
correction, subtraction of bias and dark current frames, 
linearity correction, correction for readout artifacts (the 
"bars" anomaly) , division by an appropriate flat field im- 
age, photometric calibration, cosmic ray identification, 
and combination of MULTIACCUM frames into a single 
calibrated image. No correction was made for the tem- 
perature dependence of the dark and flat-field reference 
files. These basic calibrated images (two images per filtcr 
per object) were used for the analysis. 

Figure [3] displays subsections of the calibrated imag- 
ing data for each binary in the five filter bands. Two 
components are resolved in each image, although the 
increase in PSF size in the longer wavelength data re- 
sults in better PSF separation in the F110M and F110W 
images as compared to the F170M images. There are 
notable differences in the relative component fluxes be- 
tween these two sources. The eastern component of 
2MASS J0850+1057 is significantly fainter than the west- 
ern component in all five filter bands, while the two 
components of 2MASS J1728+3948 are roughly equal in 
brightness. We measured aperture photometry for the 
combined light of each source using the IRAFn PHOT 
routine, employing a 15 pixel radius aperture centered 
on the brightest component and a 20-25 pixel annulus to 
measure the median background. Count rates were con- 
verted to Vega magnitudes on the CIT system using pho- 
tometric conversion parameters listed in the NICMOS 
Data Handbook version 8.C0, but without applying aper- 
ture corrections. Values are listed in Table [H where the 
uncertainties include shot noise and background uncer- 
tainties, as well as a standard 5% calibration uncertainty. 
Relative photometric measurements are discussed below. 

3. HST/NICMOS PSF FITTING ANALYSIS 

3 Image Reduction and Analysis Facility (IRAF; Tod"vl [19861 ) 
is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 

4 http : //www. stsci . edu/hst/nicmos/documents/handbooks/DataHandbookv8 
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TABLE 3 

HST/NICMOS Photometry and Astrometry 







2MASS J0850+1057 






2MASS J1728+3948 




Parameter 


AB 


A 


A 


B 


AB 


A 


A 


B 


Epoch 




2003 Nov 9 






2003 Sep 7 




F110W 


17.38±0.07 


1.15±0.03 


17.70±0.07 


18.99±0.08 


16.76±0.05 


0.322±0.015 


17.36±0.05 


17.68±0.05 


F110M 


17.49±0.06 


1.176±0.015 


17.81±0.06 


18.99±0.08 


16.81±0.07 


0.25±0.03 


17.44±0.07 


17.69±0.07 


F145M 


16.48±0.05 


1.000±0.015 


16.84±0.05 


17.84±0.05 


15.77±0.05 


0.517±0.014 


16.29±0.05 


16.81±0.05 


F160W 


15.70±0.05 


0.90±0.05 


16.09±0.05 


16.99±0.06 


15.08±0.05 


0.461±0.016 


15.63±0.05 


16.09±0.05 


F170M 


15.36±0.05 


0.89±0.07 


15.76±0.05 


16.65±0.07 


14.77±0.05 


0.462±0.014 


15.32±0.05 


15.78±0.05 


F110W-F160W 


1.68±0.09 




1.61±0.09 


1.86±0.10 


1.68±0.07 




1.74±0.07 


1.60±0.07 


F110W-F170M 


2.02±0.09 




1.95±0.09 


2.21±0.10 


1.99±0.07 




2.05±0.07 


1.91±0.07 


p (mas) 




132±5 






158±5 




(AU) 




5.0±0.8 






3.8±0.3 




PA (°) 




126.1±0.7 






68.0±0.5 





NOTE. — Photometry in Vega magnitudes on the CIT system based on conversions given in the NICMOS Data Handbook version 8.0. 



3.1. Method 

Relative photometry for the components of each binary 
were determined u sing a PSF- fittin g algor ithm similar to 
that described in iBurgasser et al.l (|2006cD . Models for 
each binary image were generated using PS Fs calculate d 
with the Tiny Tim package, version 6.3 (jKristl I1995I ). 
PSFs were generated for each filter passband (assuming 
post-cryocooler aberrations for NIC1), and at two point- 
ing positions to account for mirror zonal errors. The 
SpeX spectra of each source were used as templates for 
calculating filter passband effects. Individual PSF mod- 
els were generated for a subimage size of 3" x 3", which 
was resampled at ten times the original resolution to al- 
low for subpixel offsets. 

A binary PSF model for each NICMOS image was de- 
termined using an iterative image fitting algorithm that 
finds optimal positions and relative fluxes. Initial guesses 
were found using a simple peak detection algorithm and 
single-PSF subtraction on a 3" x 3" subframe of the 
data image centered on the binary. We then varied the 
primary position, secondary position, primary flux and 
secondary flux, in that order, to minimize the statistic 

52 = E E w » - M ij -( D - M hf (!) 

i 3 

where D is the data subframe, M is the model image, W 
is a masking image used to exclude bad pixels (W = for 
bad pixels), (D — M) is the mean difference (to account 
for residual background flux in the data) and the sum is 
performed over all pixels. The iteration step in position 
was set at 0.1 pixels (4.3 mas) in accord with the subsam- 
pling of the PSF model; the iteration in flux values was 
1%. Iterations were performed in a hierarchical recur- 
sive loop, and ceased when the fractional decrease in S 2 
was less than 10 -8 for all parameters. For each source, 
we adopted as our final relative flux values the mean 
of the two pointing frames in each filter band, and un- 
certainties include the difference in individual measure- 
ments and a 1% flux sampling uncertainty. Final sepa- 
rations and position angles^ were adopted as the mean 

5 The position angle is measured east of north, assuming a vector 
that points from primary to secondary. The HST roll angles at 
the time of observation a re included in the values, assumed to be 
accurate to within 0?003 HBouv et alj|200ll) . 



of all frames, with uncertainties incorporating measure- 
ment scatter and a 0.1 pixel sampling uncertainty. 

The close separation and overlapping PSFs of the 
components of these two systems initially raised con- 
cerns that fitting biases for the longer-wavelength im- 
ages (e.g., F170M) could lead to skewed component pho- 
tometry and colors. This was of particular concern for 
2MASS J0850+1057, whose components have a larger 
difference in brightness. To test our fitting program, we 
performed an identical analysis with simulated binary 
images generated from the model PSFs. The simulated 
images were constructed to have the same component po- 
sitions and relative magnitudes as determined from the 
PSF fitting (see below) as well as relative magnitudes in- 
ferred from unconstrained spectral fits (see Section 4.1). 
In both cases, the fitting algorithm reproduced the input 
relative magnitudes in all filter bands to within 0.01 mag, 
even for the F170M images of 2MASS J0850+1057. 

3.2. Results 

Table [3] summarizes the results of the PSF fitting anal- 
ysis and all of the NICMOS photometry. We find sepa- 
rations and position angles that are consisten t to within 
1.5a of those measured by IBouv et al.l (|2008l ) using the 
same NICMOS data. The relative magnitudes are all 
positive, indicating that the primary components are 
brighter than the secondary components at these wave- 
lengths for both sources. Relative component bright- 
nesses show distinct trends. For 2MASS J0850+1057, 
relative magnitudes decrease toward longer wavelengths, 
from AF110M = 1.176±0.015 to AF170M = 0.89±0.07, 
indicating a secondary that is redder than the primary. 
This is confirmed by the redder broadband color of the 
secondary, F110W-F160W = 1.61±0.09 and 1.86±0.10 
for 2MASS J0850+1057A and 2MASS J0850+1057B, 
respectively. In contrast, the relative magnitudes of 
2MASS ji728+3948 generally increase toward longer 
wavelengths (with the exception of the F145M band) in- 
dicating a bluer secondary. This trend is also seen in 
resolved JHK V photom etry for 2MASS J1728+3948 by 
IKonopackv et all (|2010D . 

Figure |3] compares the colors of the 
2MASS J0850+1057 and 2MASS J1728+3948 compo- 
nents to spectrophotometric estimates calculated from 
L2-T6 SpeX spectral templates (Section 4.1). In gen- 



6 



Burgasser, Bardalez Gagliuffi & Gizis 













* 


ft 






F110M 




F110M 






t 


F110W 




F110W 






* 


F145M 




F145M 






•* 


F160W 




F160W 


* 




* 


F170M 




F170M 



Fig. 2.— HST/NICMOS images of 2MASS J0850+1057 (left) and 2MASS J1728+3948 (right) in the F110M, F110W, F145M, F160W 
and F170M filters (top to bottom). Each image is 2"15 (50 NIC1 pixels) on a side with orientations on the sky indicated by the arrows 
(head points north, stem points east). 
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Fig. 3. — Segregation of L and T dwarfs with NICMOS photometry, based on spcctrophotometric colors computed from SpeX templates. 
The first three panels show trends in narrow-band F110M-F145M and F110M-F170M colors and broadband F110W-F160W color with 
near-infrared (SpeX-based) spectral type. The last two panels compare color pairs F110M-F145M versus F110M-F170M and F110W-F170M 
versus F110W-F160W. Measured colors for the individual components of 2MASS J0850+1057 and 2MASS J1728+3948 are indicated in 
red and blue, respectively, with the primaries (secondaries) indicated by solid (open) circles. In the first three panels, the components are 
shown with at their inferred near-infrared spectral types based on the spectral fitting analysis described in Section 4. 



eral, the components span the full range of colors seen in 
the L dwarf templates, although 2MASS J0850+1057B 
and (to lesser degree) 2MASS J1728+3948A are notably 
redder in both broadband and narrow-band colors. 
F110M-F145M and F110M-F170M colors, sampling 
1.4 /im H2O and 1.6 fim CH4 features, indicate that 
none of the components have spectral types TO or later; 
i.e., they are all L dwarfs. NICMOS colors alone are 
unable to distinguish L subtypes for the individual 



components. They are, however, cleanly segregated in 
F110W-F170M and F110W-F160W colors, which are 
correlated for the L and T dwarf templates. The binary 
components follow a sequence of 2MASS J1728+3948B, 
2MASS J0850+1057A, 2MASS J1728+3948A and 
2MASS J0850+1057B in this color-color plot, in order 
of increasing F110W-F170M and F110W-F160W colors. 

4. SPECTRAL TEMPLATE FITTING ANALYSIS 
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4.1. Method 

To infer the component near-infrared spectral types of 
2MASS J0850+1057 and 2MASS J1728+3948, we com- 
pared their combined-light SpeX spectra to a suite of 
empirical binary templates scaled to the observed rel- 
ative photometr y. We followed a proce dure similar to 
that described in lBurgasser e t al. (2006c), drawing from 
a uniform sample of low-resolution, high signal-to-noise 
SpeX prism spectra of late-type M, L and T dwarfs0 
The template sample was restricted to dwarfs later than 
M7 based on published optical classific ations for M and L 
dwarfs (tie d largely to the schemes of iKirkpatric k et al.l 
H991H991 and near-infrared classifications for T dwarfs 
(tied largely to the scheme of Burgasser eTaD 12006b") . 
About ~15% of the M and L dwarfs in our sample 
have only near-infrared types in the literature, based on 
vario us schemes (e.g., iReid et all l2001at iGeballe et al.l 
2002). We also computed alternate near-infrared spec- 
tral classifications dire ctly from the SpeX data using the 
IReid et al.1 (|2001af ) and lBurgasserl (|2007al ) index-spectral 
ty pe relations, following t he iterative procedure detailed 
in IBurgasser et al.l (|2010D ■ Known (resolved) binaries 
and sources noted to have peculiar spectra or highly un- 
certain classifications (a > ±2 subtypes) were purged 
from the template sample. The resulting 462 spectra 
of 438 sources were interpolated onto a common wave- 
length scale, and synthetic Vega magnitudes in the five 
NICMOS filters and MKCQ JHKK p filters were com- 
puted by convolving each spectrum and a Kurucz model 
of V ega with the appropr iate filter transmission functions 
sec ICushing et aLl^005l) . 

From these individual templates, we produced sepa- 
rate libraries of binary templates for 2MASS J0850+1057 
and 2MASS J1728+3948 by adding together appropri- 
ately scaled pairs. The pairs were initially selected to 
have secondary component types that were no more 
than 2 subtypes earlier than the primary type. This 
produced 137,212 unique combinations. We then com- 
puted an uncertainty-weighted mean relative flux scal- 
ing between the two components of each system using 
all fi ve NICMOS measurem ents and Kp photometrjjfl 
from iKonopackv et all (|2010D ; see Table [lj After scaling 
the spectra, we required that relative synthetic magni- 
tudes in all six filter bands agree with measured values 
to within 3a. This substantially reduced the number of 
"allowable" templates to 19,042 for 2MASS J0850+1057 
and 1389 for 2MASS J1728+3948. 

Comparisons were then made between the source 
spectral data and allowed binary templ ates using the 
weighted chi-squared statistic defined in ICushing et al.l 



Ull), 



6 These da ta were c o mpiled fromlBurgasser e t al.l j2D03.l200 6aHH 
|2007|. 12008a b, 2010); Cruz et all (20 04); Burgasser & M cElwainl 
(20061): IChiu et ; al.l 112009); IMcElwai n fc Burgasserl 120061)1* 
IReid et all (2009); iBurgasserl (2007a. b e) ; ILiebert fc Burgasser] 
(20071); ILooper et aTj (20071. I20tl8bl); ILuhman et all (20071 ); 
Siegler et all H2007I) ; ISheppard fe Uushinal (2009) and 
""chrm"aT"e*t"*a l (2010a), and are available at the SpeX Prism Spec- 



tral Libraries website, http://www.browndwarfs.org/spexprism 

7 M auna Kea Observatory filter s ystem; see Tokunag a et all 
(2002) and Simons & Tokunaga (2002). 

8 Konopacky ct al. (2010) also report resolved J- and H- 
band photometry for 2MASS J1728+3948, AJ = 0.32±0.02 and 
AH = 0.45±0.02. These measurements are redundant — and in 
agreement— with broad-band F110W and F160W NICMOS pho- 
tometry, and were therefore not used in the analysis. 



{A} 



,[A] 



D[X] - aT[X] 



o-d[X] 



(2) 



Here, D[X] and T[A] are the data and template spectra, 
respectively; o\d[A] is the uncertainty spectrum of the 
data; w[X] is a vector of weights satisfying X){a} w l^\ = 1> 
a is a scaling factor that minimizes Gk (see Equation 2 in 
ICushing et al.ll2008D ; and the sum is performed over the 
wavelength ranges {A} = 0.95-1.35 /im, 1.45-1.8 /xm and 
2.0-2.35 /im to avoid regions of strong telluric absorption. 
We adopted th e same weig hting scheme as that used in 
ICushing et all (|2008l ) and (Burgasser et al.l (|2010D . with 
each pixel weighted by its spectral width (i.e., Wi = AA^). 

Templates which provided minimum values of Gk were 
deemed to be the best fits. However, it was gen- 
erally the case that other templates gave Gk values 
only slightly larger than this best-fit template, and are 
therefore statistically equivalent. We therefore deter- 
mined average component parameters (spectral types 
an d relative magnitudes) and their uncertainties follow- 
ing (BurgasseFeT^y (|2010l Equations 4-6), with an effec- 
tive degrees of freedom v = 253 used for weighting the 
parameters with the F-test statistic. 

4.2. Results 

The resulting best-fitting binary templates are listed 
in Table 21 while Figures [1] and [5] display the four best 
fits to 2MASS J0850+1057 and 2MASS J1728+3948, re- 
spectively. These templates reproduce the observed data 
rather well, including the overall spectral slopes, shapes 
of the JHK flux peaks, and depths of the H2O, CO and 
FeH absorption bands. By design, each of the spectral 
combinations listed also reproduce the relative magni- 
tude measurements to within the observational uncer- 
tainties. Notably, the binary template fits are statis- 
tically superior to equivalent single template compar- 
isons. Best-fit Gk values for these single templates were 
2.20 for 2MASS J0850+1057 ( comparing t o the L5/L6.33 
2MASSW J1326201-272937; iGizisI [2002) and 3.90 for 
2MASS J1728+3948 (comparing to the L5:/L5 2MASS 
J06244595-4521548; [Reid_elalj|2008) , differing at the > 
99.9% significance level based on the F-test statistic. The 
significantly improved fits by the binary templates indi- 
cate that they are (as expected) better representations 
of the combined light spectra of these sources. 

The spectra of the best-fit template components also 
replicate photometric trends from HST and LGSAO 
imaging. For 2MASS J0850+1057, each of the tem- 
plates in Figure [4] has a secondary that is distinctly red- 
der than its primary. These secondaries — 2MASSW 
J1553214+210907 (L5.5 /L6.1:, J - K s = 2.03±0.19: 
IKirkpatrick eTaTl fl999l) . 2MASSI J0028394+150141 
(L4.5 /L5.2::, J — K s = 1.95±0.13; IKirkpatrick "eTall 
l2000T ) and SDS Sp J010752.33 +004156.1 (L8/L6.0, J-K s 
= 2.12±0.07; IGeballe et all l2002f ) are among the red- 
dest L dwarfs currently known. In contrast, the pri- 
maries of these templates have J — K s colors that are 
in line with median values for their spectral types (e.g., 

9 Hereafter, we list both literature and SpeX spectral types for 
each template; uncertain types are indicated with ":" for ±1 sub- 



type, 



for ±2 subtypes. 
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TABLE 4 

Best Fitting Binary Templates 



Primary 


SpT 


SpT 


Secondary 


SpT 


SpT 


G k 


Relative 




(Lit.) a 


(SpeX) b 




(Lit.) a 


(SpeX) b 




Weight 


2MASS J08503593+1057156 


SDSS J151506.11+443648.3 


L7.5±1.5 d 


L6.9 


2MASSW J1553214+210907 


L5.5 


L6.1: 


1.46 


1.00 


SDSS J151506. 11+443648.3 


L7.5±1.5 d 


L6.9 


2MASSI J0028394+150141 


L4.5 


L5.2:: 


1.57 


0.28 


SDSS J151506.11+443648.3 


L7.5±1.5 d 


L6.9 


SDSSp J010752.33+004156.1 


L8 


L6.0 


1.62 


0.20 


SDSS J171714.10+652622.2 


L4 


L5.8: 


SDSSp J010752.33+004156.1 


L8 


L6.0 


1.85 


0.03 


SDSS J171714.10+652622.2 


L4 


L5.8: 


SDSS J080959. 01+443422.2 


L6 d 


L6.2: 


1.88 


0.02 


SDSS J171714.10+652622.2 


L4 


L5.8: 


SDSSp J132629.82-003831.5 


L8: 


L6.1: 


1.94 


0.01 


2MASS J03101401-2756452 


L5 


L6.0: 


SDSSp J010752.33+004156.1 


L8 


L6.0 


1.95 


0.01 


Weighted Mean 


L7.5±0.0 


L6.9±0.0 


Weighted Mean 


L5.7+1.2 


L5.8+0.4 






2JV1ASS J 17281150+3948593 


2MASS J01443536-0716142 


L5 


L3.9 


SDSS J104409. 43+042937.6 


L7 a 


L7.1 


2.41 


1.00 


2MASSW J2224438-015852 


L4.5 


L4.0: 


SDSS J104409. 43+042937.6 


L7 d 


L7.1 


2.58 


0.29 


2MASS J03185403-3421292 


L7 


L6.5 


2MASS J09054654+5623117 


L5 


L5.6: 


2.64 


0.23 


2MASS J01443536-0716142 


L5 


L3.9 


2MASS J23254530+4251488 


L8 


L7.1 


2.75 


0.15 


2MASS J03185403-3421292 


L7 


L6.5 


2MASSI J1305410+204639 


L4: 


L6.0: 


2.92 


0.06 


2MASS J01443536-0716142 


L5 


L3.9 


SDSSp J003259. 36+141036.6 


L8 d 


L7.9: 


3.12 


0.02 


Weighted Mean 


L5.3±1.0 


L4.6+1.1 


Weighted Mean 


L6.5+1.1 


L6. 8+0.6 







a Optical spectral type from the literature unless otherwise noted. 

b Classifications base d on H2O and CH4 in dex/spectral type relations defined in Burgasscr (2007a), following the iterative 
procedure detailed in Burgasscr ct al. (2010). 

c Statistical weight assigned to template parameters in calculation of means and uncertainties. This weight corresponds to the 
F-test probability dis tribution function for the ratio Gfc/min(Gfc) and effective degrees of freedom v = 211 (see Equation 4 in 
Burgasscr ct al. 2010). 

d Near-infrared classification from the literature. 

J104409.43+042937.6 (X7/L6.9: iKnatro et al.ll27)0l and 
2MASS J23254530+4251488 rL8/L7.3:: ICruz et al.ll2007T) 
show weak signatures of 1.6 and 2.2 \ra\ CH 4 absorption 
in their near-infrared spectra. These features suggest 
that 2MASS J1728+3948B is at the threshold of the L 
dwarf/T dwarf transition, despite its mid-L near-infrared 
spectral type. We adopt mean component types of L5 
and L6.5 for this system. 

4.3. Assessment of Systematic Effects 

To assess the robustness of our results, we conducted 
the same template fits for three subsets of relative pho- 
tometry: the NICMOS data alone, and single filter scal- 
ings with F110W and K p photometry. The NICMOS 
fits produce little change in the inferred spectral com- 
ponents, with 2MASS J0850+1057 still hosting a later- 
type primary (L7 + L6 components based on both liter- 
ature and SpeX classifications) and 2MASS J 1728+3948 
having a primary that is typed earlier than the com- 
bined light spectrum (L5 + L6.5 for literature classi- 
fications, L4 + L7 for SpeX classifications). However, 
the mean relative K p magnitudes from these fits — 
0.950±0.015 for 2MASS J0850+1057 and 0.54±0.05 for 
2MASS J1728+3 9 48 — differ from the measurements of 
IKonopackv et al.l (|2010f ) at the la and 1.5a levels, re- 
spectively. Inferred brightnesses diverge even more dra- 
matically when a single filter is used to scale the tem- 
plates. FllOW-scaled templates for 2MASS J0850+1057 
indicate component types of L6 + L9:, with a signifi- 
cantly reduced best-fit Gk value (1.12 versus 1.46, dis- 
tinct at the 98% confidence level). These types ar e 
more in line with prior estimates (e.g. JReid et al.lfeOOlbl) . 
but the inferred component magnitudes differ signifi- 
cantly from measured values: AF160W disagrees at the 
8a level, while AK p disagrees at the 6.5a level. K p - 



ISchmidt et aHl2010bD . What is more surprising, how- 
ever, is that the best-fit secondaries are generally of ear- 
lier type than the best-fit primaries. This is true for both 
literature and SpeX classifications. The mean compo- 
nent classifications reflect this: L7.5 and L5.5±1.2 based 
on literature classifications and L7 and L6 based on SpeX 
classifications; we adopt the latter combination as the 
uncertainty- weighted means. This apparent reversal in 
classifications is particularly remarkable given the 0.8- 
1.2 magnitude brightness difference between the primary 
and secondary. There is no evidence that either compo- 
nent is a T dwarf. 

For 2MASS J1728+3948, the best-fit secondaries are 
consistently bluer than the primaries, again in line with 
photometric trends. Notably, the components of the bi- 
naries shown in Figure [5] have comparable fluxes in the 
0.9-1.0 Atm region, supportin g the F1042M flux reversal 
reported bv lGizis et al.l (|2003[ ) , and diverge toward longer 
wavelengths. For these templates, it is the primary 
components that are unusually red: 2MA SS J01443536- 
0716142 (L5/L4.2, J - K s = 1.92±0.03; iLiebert et all 
12001 . 2MASSW J2224438-015852 (L4. 5/L4.5:, J - 
K s = 2.05±0.04; IKirkpatrick eHd! [2000h and 2MASS 
J03185403-342129 2 (L7 /L6.4, J — K s = 2.06±0.07; 
IKirkpatrick et all 120081 ). The mean spectral type of 
2MASS J1728+3948A, L5.5±1.0 from literature clas- 
sifications and L4.5±l.l from SpeX classifications, is 
also remarkable for being considerably earlier than the 
combined-light L7 optical classification. It is typically 
the earlier-type component that dominates optical flux 
in L and T dwarf binaries. In this case, it appears that 
the comparable brightnesses of the two components at 
red optical wavelengths produces a "blended" combined- 
light spectral type. The secondaries of the best-fit 
templates have normal colors, although both SDSS 
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G k = 1.46 

SDSS J151506.11+443648.3 L7.5+/-1.5 (L6.9) 
2MASSW J1553214+210907 L5.5 (L6.1:) 




G k = 1.57 
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Fig. 4.— The four best template fits for 2MASS J0850+1057 as constrained by HST/NICMOS and LGSAO K p photometry. In each 
panel, the source spectrum is indicated by the black line, the source noise spectrum by the grey line, the best-fitting template by the 
green line, the best-fitting primary by the red line and the best-fitting secondary by the blue line. Template spectra (and appropriately 
scaled primaries and secondaries) are shown at their minimum G^ scalings relative to the source spectrum, which is normalized at its peak 
spectral flux. Template source names, literature classifications, SpeX-based classifications (in parentheses) and fit values are indicated. 
Note that the L7.5±1.5 literature spectral type for SDSS J151506. 11+443648.3 is based on near-infrared data; all other literature types are 
based on optical data. 



scaled templates for 2MASS J0850+1057 produce com- 
parably large discrepancies in inferred NICMOS photom- 
etry. Single-filter template fits to 2MASS J1728+3948 
are somewhat more robust, producing similar compo- 
nent types as the baseline template sample (L5+L7 for 
F110W fits, L5+L7.5: for K p fits), although inferred rel- 
ative magnitudes still deviate by up to 2a. 

We conclude that multi-band photometry is essen- 
tial to accurately and precisely constrain the component 
properties of these binary systems, and particular care 
must be taken when drawing conclusions from single- 
filter measurements. Nevertheless, the general agreement 
in component types inferred from different multi-band 
subsets (e.g., NICMOS+if p , NICMOS only) indicates 
that our results are robust. 

5. ABSOLUTE MAGNITUDE/SPECTRAL TYPES AND 
COLOR TRENDS 

Both 2MASS J0850+1057 and 2MASS J 1728+3948 
have parallax distance measurements from IVrba et all 
(|200l . so it is possible to compare component abso- 
lute fluxes to those of comparable field dwarfs. Fig- 



ure IH1 displays combined light and component absolute 
MKO J and K magnitudes versus near-infrared spec- 
tral type and J — K color. The MKO magnitudes from 
2MASS J0850+1057 a nd 2MASS J17 28+3948 are based 
on measurements from iLeggett et al.l (|2002T ) for the for- 
mer, and synthesized from 2MASS photometry and SpeX 
spectroscopy for the latter. Component JK fluxes were 
determined from the template fits, synthesized from the 
various best-fit template spectra and combined following 
the same weighting scheme as the average spectral types 
(Table [5]). We use the near- infrared spectral types calcu- 
lated in Section 2.2. Th e comparison sample was drawn 
from the compilation of ILeggett et al.l ()2010[ ). where we 
have either used published near-infrared spectral types 
listed in that study or types calculated from SpeX spec- 
tra (where available) following the same spectral index 
method described above. 
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G k = 2.41 

2MASS J01443536-0716142L5 (L3.9) 
SDSS J104409.43+042937.6 L7 (L7.1) 




Wavelength (Jim) 



G k = 2.58 

2MASSW J2224438-015852 L4.5 (L4.0:) 
SDSS J104409.43+042937.6 L7 (L7.1) 




1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 

Wavelength (\lm) 



G k = 2.64 

2MASS J03185403-3421292 L7 (L6.5) 
2MASS J09054654+5623117L5 (L5.6:) 




Wavelength (|lm) 



G k = 2.75 

2MASS J01443536-0716142 L5 (L3.9) 
2MASS J23254530+4251488 L8 (L7.1) 




Wavelength (flm) 



Fig. 5. — The four best template fits for 2MASS J1728+3948, displayed as in Figure |4] Note that the L7 literature spectral type for 
SDSS J104409. 43+042937. 6 is based on near-infrared data; all other literature types are based on optical data. 
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Fig. 6.— (Top panels) Absolute MKO J and K magnitudes versus spectral type for the 2MASS J0850+1057 (red) and 2MASS J1728+3948 
(blue) components. These are compared to single (open symb ols) and binary L and T dwarfs (filled symbols) with absolute magnitude 
uncertainties < 0.3 mag, as compiled by Lcggctt ct al. (2010). All sources are plotted according to their near-infrared classifications, 
either pu b lished or computed from SpeX spectroscopy. Also shown are the absolute magnitude/spectral type relations quantified in 
Liu ct al. (2006, dashed lines); both "bright" and "faint" relations are shown. Combined light photometry for 2MASS J0850+1057 and 
2MASS J1728+3948 are connected to component values. (Bottom panels) Absolute magnitudes versus MKO J — K color based on the 
same sample. 
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Focusing first on the absolute magnitude/spectral type 
comparisons, we find that three of the components — 
2MASS J0850+1057B, 2MASS J1728+3948A and 
2MASS J1728+3948B— roughly cluster, with ab- 
solute magnitudes similar to other L6-L8 sources. 
2MASS J1728+3948A, and to a lesser degree 
2MASS J1728+3948B, are somewhat underlumi- 
nous for their types, particularly at J-band where 
they fall ps0.5 mag belo w the absolute magnitude 
relations of lLiu et all (p006). In contrast, the primary of 
2MASS J0850+1057is «1 mag brighter than comparable 
L6-L8 dwarfs, or equivalently has a classification 3 sub- 
types too late for its measured absolute brightness. This 
component is even marginally brighter than combined 
light photometry for 2MASS J1728+3948. As such, 
the discrepancy between brightness and spectral type 
between the components of the 2MASS J0850+1057 
system appears to be rooted to the unusually late 
spectral type and/or overluminosity of its primary. 

In the near-infrared color-magnitude diagrams 
(CMDs), all four components follow the locus of field 
L dwarfs, tracing the inflection in J — K color and Mj 
magnitude at the end of the L sequence. In particular, 
2MASS J0850+1057B and 2MASS J1728+3948B bridge 
the turnover to bluer near-infrared colors, and the latter 
appears to be the faintest L dwarf with a parallax 
measurement currently known. Its location on the 
near-infrared CMDs supports our hypothesis that it 
is on the cusp of becoming a T dwarf (Section 4.2). 
The components of 2MASS J0850+1057 are no longer 
outliers in these plots, in contrast to prior results (e.g., 
IVrba et alJl200l . 2MASS J0850+1057A in particular 
does not stand out as unusually bright, due largely 
to the near-vertical locus of L dwarfs in near-infrared 
CMDs. 

For completeness, we also determined luminosity, T e ff 
and mass estimates for the components of these bi- 
naries using their inferred absolute magnitudes and 
near-infrared spectral types. Bolometric luminosities 
were calculated using if-band bolometric corrections 
(BCk) derived from the BCk /spectral type relation of 
iDupuy fc Liu! ()2009f ). Note that these values may be 
systematically too high for the unusually red sources 
2MASS J0850+1057B and 2MASS J1 728+3948A. T P . ff s 
and m asses were estimated using the lSaumon fc Marlevl 
(2008) cloudy evolutionary models, based on the de- 
rived luminosities and age estimates listed in Ta- 
ble [U for 2MASS J1728+3948, we assumed an age 
range of 1.5-7.5 Gyr, the upper limit based on 
the small Vt an of this source. The inferred T e ffS 
for 2MASS J1728+3948A, 2MASS J0850+1057B and 
2MASS J1728+3948B are consistent with pri or esti- 
mates for L5-L6.5 dwarfs (Stephens et al. 2009), while 
the brightness of 2MASS J0850+1057A makes it «300 K 
warmer than comparable L7 dwarfs. We discuss this 
component in detail in Section 6.2. 

6. DISCUSSION 

6.1. Clouds and Classification 

Our photometric and spectroscopic analyses 
of the components of 2MASS J0850+1057 and 
2MASS J1728+3948 have revealed several unusual 
traits, particularly in their primaries: under- and overlu- 



minous fluxes and unusually early and late near-infrared 
spectral types. These peculiar traits can be related to 
their unique atmospheric properties. 

In the case of 2MASS J1728+3948A, we hypothe- 
size that condensate cloud effects are responsible for 
shifting this component toward both an earlier near- 
infrared spectral classification and toward slightly fainter 
J-band fluxes. Grey extinction from condensate cloud 
grains in L dwarf photospheres dominate the opacity 
at the J- and ii-band flux peaks, as these windows 
in molecular gas opacity probe deeper into the atmo- 
sphere and sample a larger column depth of cloud ma - 
terial ([Ackerman fc Marlevl 120011 : iBurrows etall [2006) . 
The K-b&nd peak, on the other hand, is modulated by 
both cloud op acity and collision-induced Hi absorption 
(|Linsky|[l969f ). As a result, greater condensate opacity 
tends t o produce redder J — K colors and fainter J-band 
fluxe s ([Allard et al.ll2001tlMarley et al.ll200l iTsuji etall 
l2004f ) . In addition, contrast in molecular absorption 
bands is reduced, particularly for the near-infrared H 2 
and FeH ba nds to which near-infrared schemes ar e com- 
monly tied ([Reid et al.ll200Tat iGeballe et alll2002l ). Veil- 
ing of these features c an skew near-infrared classification s 
toward earlier types ([Stephens 20031 : iKnapp et al.ll2004T ) . 
Hence, thicker condensate clouds leads to redder near- 
infrared colors, reduced flux at J-band, and systemati- 
cally earlier near-infrared classifications. 

These conditions accurately reflect the properties 
of 2MASS J1728+3948A, and to a lesser degree 
2MASS J0850+1057B. The best-fit primary compo- 
nents for 2MASS J1728+3948 are consistently red, mid- 
type L dwarfs whose near-infrared spectral types are 
consistently earlier than their optical types (Table [4]). 
These include the templates 2MASS J01443536-0716142 
(hereafter 2MASS J0144-071 6; L5 optica l clas sifica- 
tion, L4 SpeX classification; iLiebert et all 120031 ) and 
2MASSW J2224438-015852 (hereafter 2MASS J2224- 
0158; L4. 5 optical classification, L3 . 5 near-infrared clas - 
sification; iKirkpatrick et al.l 120001: IKnapp et all 120041) . 
Both of these sources exhibit indications of thick clouds, 
based on the detection of linear polarization in the 
case of 2MASS J0144-071( g] (0.58±0.19% at /-band; 
iZapatero Osorio et al.ll2005l ). and spectral model fits and 
pronounced silicate grain ab sorption at 9—11 jjm in the 
case of 2MASS J222 4-0158 (ICushing et all 120061 120081: 
IStephens et all 120091 ). Unusually red dwarfs such as 
these have bee n found to have cool er T e ffS for their 
spectral types ([Stephens et all 120091 ) . which can con- 
tribute to fainter magnitudes. The older age of the 
2MASS J1728+3948 system, based on the absence of Li I 
absorption, argues that thick condensate clouds, rather 
than l ow surface gravity, gives rise to its un usually red 
color ([Allers et alJl2007t Hooper et al.ll2008bl ). 

Thick clouds in 2MASS J 1728+3948 A may also ex- 
plain its comparable brightness at 1 /xm compared 
to its later-type companion. As discussed above, 
2MASS J1728+3948B appears to be at the threshold of 
the L dwarf/T dwarf transition, a phase in which con- 
densate clouds are inexplicably dispersed and J-band 
fluxes increase. In a sa mple of unresolv e d L/T tran- 
sition binary candidates, iBurgasser et all (|2010l) found 

10 2MASS JQ144 -0716 was also detected in an optical flare by 
ILiebert et all l(20ul ). 
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TABLE 5 

Adopted Component Parameters for 2MASS J08503593+1057156 and 
2MASS J17281150+3948593 



Parameter 



2MASS J0850+1057 
A AA a B 



2MASS J1728+3948 
A B 



NIR Spectral Type 


L7 


L7+L7 


L6 


L5 


L6.5 


Mj 


13.62±0.35 


14.37±0.35 


14.75±0.35 


14.60±0.19 


14.90±0.19 


M K 


11.82±0.35 


12.57±0.35 


12.77±0.35 


12.48±0.18 


13.07±0.18 


J - K 


1.80±0.04 




1.98±0.05 


2.12±0.09 


1.84±0.09 


logL 6oi /L b 

Te// (K) c 


-4.13±0.14 


-4.43±0.14 


-4.52±0.14 


-4.42±0.08 


-4.63±0.08 


1770±200 


1480±180 


1400±170 


1570±80 


1340±70 


Mass (M Q ) C 


0.06±0.02 


0.10±0.03 


0.05±0.02 


0.075±0.007 


0.066±0.008 



a Parameters assuming two equal-mass components. Mass is given for combined pair. 

b Luminosity based on Mfc magnitudes and X-band bolometric corrections (BCpc) from 

IDupuv fc Li u (2009). Stated ler errors include ±0.5 uncertainty in spectral type, absolute 

magnitude uncertainties, and 0.08 dex scatter in the Dupuy & Liu BC/f/spectral type 

relation. 

c T e ff and mass estimates based on inferred luminosities and ages listed in Table [T] 
(2MASS J1728+3948 is assumed t o be 1.5-8 Gyr), combined with cloudy evolutionary mod- 
els from Saumon & Marlcy (2008). Stated la errors include uncertainties in luminosities 
and spread in age estimates. 



that sources with comparable component types but red- 
der primaries showed a more pronounced flux reversal 
than sources with normal or blue primaries. They argued 
that J-band fluxes in the primaries of these systems were 
more suppressed. We may be seeing a similar effect in 
the 2MASS J1728+3948AB pair. Alternatively, we may 
be observing the tops of thick clouds that are constrained 
to the same temperatur e layer in both sources (i.e., T cr 
— constant; iTsuiil I2005T ). These possibilities should be 
explored with detailed modeling of resolved component 
spectra, rather than spectral templates. 

6.2. Is 2MASS J0850+1057 a Young Triplet? 

While thick clouds can explain the unusual faintness 
of 2MASS J1728+3948A, thin clouds do not readily ex- 
plain the unusual brightness of 2MASS J0850+1057A. 
Thin-cloud L dwarfs, also known as unusually blue L 
dwarfs (UBLs), do exhibit discrepa ncies between opti- 
cal and near-inf r ared sp ectral types (jKnapp et al.ll2003 : 
iBurgasser et al.1 1200 8b). There is also evidence that 
UBLs have systematically larger T e //s and /or absolute 
fluxe s for their spectral classifications dStephens et al.l 
I2009D . However, 2MASS J0850+1057A has a normal 
near-infrared color for its spectral type, and best-fit tem- 
plates exhibit none of the spectroscopic hallmarks of a 
UBL (e.g., deep H 2 Q and FeH bands, blue near-in frared 
SEP; IBurgasser et alj[20(58bt iSchmidt et al.ll2010al) . It is 
therefore unlikely that this source is a UBL with unusu- 
ally thin clouds. 

Youth may play a role in the unusual brightness of 
2MASS J0850+1057A, arising from the enlarged radius 
of a still-contracting brown dwarf. The radius of a 1500 K 
brown dwarf is 25% larger at an age of 250 Myr - 
the minimum bound cited by iFahertv et all (|2010l) - 
than at 2 Gyr, corresponding to an increase in bright- 
ness of roughly 0.5 mag. Given the current uncertain- 
ties in the distance of this source, such a shift may be 
sufficient to move 2MASS J0850+1057A back onto ab- 
solute magnitude/spectral type tracks. However, assum- 
ing that the two components of this system are coeval, 
2MASS J0850+1057B would also have to be overlumi- 
nous by roughly the same factor, which does not appear 



to be the case. Moreover, such a correction fails to ex- 
plain the significant brightness difference between these 
two comparably-classified sources. 

We propose an alternative explanation: 
2MASS J0850+1057A is itself an unresolved, near- 
equal mass binary. Such a scenario would explain how 
this component could appear both brighter and warmer 
than its equivalently-typed companion, but have an 
otherwise normal spectral energy distribution. The 
components of an equal-mass 2MASS J0850+1057A pair 
would have absolute magnitudes, luminosities and T e f fS 
fully consistent with empirical trends (Table [5]). The 
fact that existing high angular resolution images have 
not resolved this source requires an angular separation 
<50-100 mas, or a projected separation <2-4 pc at 
this distance of this system. In fact, long-term dynamic 
stability requires an even tighter binary. A hierarchical 
triple is generally stable if the ratio of outer periastron 
and inner apastron distances: 



satisfies 



Y > Y r ' 



Y 



3.7 

T/3 
lout 



t(l - e out ) 
t (l + e in ) 



2.2 



o 1/3 

Hout 



1 



+ 



1.4 



1/3 

^OUt 



1/3 1/3 
lout 



1 



(3) 



(4) 



> 1 and 
mass ratios. 



(Eggl eton fc Kiseleval[l995T ). Here, q m = f| 

Qout = Mi m 3 M2 are the inner and outer 
Assuming that all three components have nearly equal 
masses, and that the inner and outer orbits have eccen- 
tricities e = (0.5), Equation 2] requires a limit on the 
inner orbit semimajor axis of a,i n < 0.2a ou t (<0.07a ou t) 
and hence a.;„ < 25 mas (8 mas), or 1 A U (0.3 AU), based 
on the semimaj or axis dete r minat ion of lKonopackv et al.l 
(|2010T l and the lVrba etHI (|2004| ) parallax. 

Such a tight separation is not unusual for brown 
dwarf multiples; brown dwarf spectroscopic binaries 
with compar able separations have already been iden- 
tified (e.g., iBasri fc Martini 119991: iBlake et al.l 120081: 
Uoergens et al.H2010D . Moreover, three other tight, hi- 
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erarchical, brown dwarf triple candidates identif ied in 
the literature— Gliese 569BCD. (iSimon et al.l l2006ft. DE- 
NIS J0205-1159ABC (|Bouv et al.ll2005D . and Kelu 1ABC 
(jStumpf et al.ll2010h — exhibit evidence that one compo- 
nent is an unresolved pair, based on radial velocity vari- 
ations, PSF-fitting residuals and spectroscopic features, 
respectively. This scen ario is also consistent wit h orbital 
mass constraints from iKonopackv et all (|2010() , as the 
estimated total mass of a triple 2MASS J0850+1057 sys- 
tem, 0.15±0.04 Mq, is closer to the mean (but weakly 
constrained) value of 0.2 M Q found in that study. In- 
deed, tighter constraints on the total mass of this system 
from ongoing astrometric monitoring may affirm or re- 
fute the presence of a third body. The triple hypothesis 
can also be tested though radial velocity monitoring; a 
<5 km s _1 (<8 km s _1 ) line shift arising from a pair 
of 0.05 M Q brown dwarfs separated by 1.0 AU (0.3 AU) 
can be read ily detected with current near-infrar ed instru- 
mentation (|Blake et al.ll2007HBean et al.ll2010l) . 

If 2MASS J0850+1057A is confirmed as a binary, 
it would complete a remarkable, young, low-mass, hi- 
erarchical quintuple system with the double M dwarf 
NLTT 20346 AB, encompassing 4 orders of magnitude in 
separation and composed entirely o f objects less massive 
than 0.15 M (jFahertv et al.ll2010h . 

7. SUMMARY 

We have presented photometric and spectro- 
scopic analyses of the late-type L dwarf binaries 
2MASS J0850+1057 and 2MASS J1728+3948, aimed at 
assessing component spectral types, absolute magnitudes 
and near-infrared colors. Multi-band HST/NICMOS 
photometry have revealed distinct trends in the relative 
colors of these two systems, with 2MASS J0850+1057B 
being redder than its primary and 2MASS J1728+3948B 
being bluer. Neither secondary exhibits narrow-band 
colors consistent with being a T dwarf. These re- 
sults are borne out in spectral template fits, using 
NICMOS and K p resolved photometry, which also 
determine component near-infrared spectral types of 
L7 + L6 for 2MASS J0850+1057 and L5 + L6.5 
for 2MASS J1728+3948. The early classification of 
2MASS J1728+3948A, its relative faintness at J, and 
its unusually red color can be explained by the presence 
of thick condensate clouds in its photosphere. The 
secondary of this system, in contrast, may be losing 
its photospheric cloud deck as it transitions onto the 
T dwarf sequence. For 2MASS J0850+1057, the sur- 
prisingly later spectral type of its bright primary may 
stem from youth (inflated radius) and/or unresolved 
multiplicity. The latter hypothesis, which would make 
2MASS J0850+1057 part of a low-mass hierarchical 
quintuple, can be tested through ongoing astrometric 
monitoring and/or resolved spectroscopic monitoring to 
search for RV variations. 

As two resolved (or partly-resolved) coeval systems 
spanning the end of the L dwarf sequence and exhibiting 
a broad range of cloud properties, 2MASS J0850+1057 
and 2MASS J1728+3948 remain important laborato- 
ries for studying cloud formation and evolution in low- 
temperature atmospheres. Improved parallactic distance 
measurements — including resolution of current distance 
discrepancies for 2MASS J0850+1057 — resolved compo- 
nent spectroscopy, and ongoing photometric and astro- 



metric monitoring will aid in characterizing the clouds, 
spectral properties and multiplicity of these benchmark 
brown dwarf systems. 
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